NMR spectroscopy and computing beyond myth and fiction 
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Starting from the physical algebra, namely the algebra of operators corre- 
sponding to macroscopic measurements, a thorough classical theory of NMR 
spectroscopy is constructed. Its effectiveness in emulating quantum comput- 
ing stems from the exponential growth of the classical phase space dimension 
with the number of magnetic nuclei per molecule, this solving a standing 
scientific dispute. As a by-product, the differential equations for the evolu- 
tion of macroscopic variables give an exact setting for simulation of NMR 

spectroscopy both in general and for quantum computing emulation. 
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qh| At room temperature relaxation and decoherence times for the spin degrees of freedom of 



magnetic nuclei of suitable molecules in liquid solutions put in strong magnetic fields of some 



^ I Tesla are exceedingly long, even several seconds [|I| . In concert with the presence of chemical 
shifts, giving rise to small differences in Larmor frequencies for nuclei of the same kind, and 
of the weak J-coupling between them, this would make such systems ideal candidates for 
quantum computing [|| if only one could control these degrees of freedom. Even though 
this is not realistic for a single molecule (s.m.), it is achieved by NMR spectroscopy in a 
statistical sense for macroscopic samples. NMR experiments are traditionally analyzed in 
terms of statistical ensembles of s.m.'s and, if the J-coupling is ignored, each molecule is 
equivalently described either classically or quantum mechanically, due to the linearity of 
the equations of motion |^ . This same linearity makes the difference between the ensemble 
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viewpoint and the real physical situation, where measurements are done on the whole sample, 
partially irrelevant. On the other hand the introduction of the J-coupling in order to describe 
the so called multiple resonances, within the ensemble approach seemingly forces to use the 
quantum description 

This description in terms of an ensemble of quantum systems is taken as the starting 
point of the nice idea of emulating quantum computing by NMR spectroscopy [|],^. In 
particular this is done using the fact that the traceless part of a single particle thermal state 
can be identified with the traceless part of a pure state, which gives rise to the notion of a 
pseudo-pure state. 

The above setting is physically unsatisfactory since the macroscopic high temperature 
character of the states involved and the classical nature of the measurements performed, 
namely the absence of wave function collapse, are an unambiguous indication that the 
soundest description is classical. On the other hand the quantum ensemble description, 
in addition to its lack of physical cogency in treating the coherent electromagnetic field pro- 
duced in modern pulse NMR interferometry, is a source of considerable confusion. In fact, 
as reported in a recent paper [Q], the quantum-mechanical nature of processes involving such 
macroscopic pseudo-pure states is a matter of debate. In paper |^, for instance, it is argued 
that, although since the first proposals to use NMR for quantum computation, there has been 
surprise about the apparent ability to perform quantum computations in room-temperature 
thermal ensembles, entanglement analysis suggests that current NMR experiments should be 
considered as simulations of quantum computations, rather than true quantum computation. 
In other papers on the contrary one can find explicit references to the quantumness of NMR 
{Perhaps the most interesting aspect of this experiment is that the experimental results can- 
not be explained by a classical model of two interacting spins ^) and to the presence of 
entangled states {Because the N spins in each molecule may be in entangled quantum super- 
position states, our computer is a quantum one Then, while the single particle quantum 
ensemble paradigm turns out to be a fruitful viewpoint in analyzing NMR spectra for struc- 
ture analysis in chemistry, a more transparent physical setting is needed if one wants to 
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address more fundamental questions. This is witnessed by recent attempts to analyze the 
structure of the space of density matrices corresponding to the aforementioned pseudo-pure 
states p|, p!3|JT^ , to construct classical probabilistic models to mimic NMR experiments 
and to analyze models for ensemble computation in general |I0[| . 

The aim of this paper is to show that, by sticking to the physical algebra [|IT|, namely 



the operator algebra physically relevant to the system in consideration, a somehow definitive 
answer to the fundamental questions raised by the NMR realization of quantum algorithms 
can be given. In fact it turns out that the correct collective quantum description, like in the 
case without J-coupling, is equivalent to the corresponding classical one, where expectation 
values replace operators. The physical description of the whole sample can be recast in 
such a way to be formally equivalent to the one previously used for a s.m. and allowing 
quantum computing emulation. In particular the thermal state of the physical algebra, 
which is a convex combination of states belonging to inequivalent irreducible representations 
T5| , plays the role of a pseudo-pure state. As a relevant by-product, the present approach 



leads, without any approximation, to a system of linear differential equations for classical 
observables and classical inner variables, which may be of use in the analysis and in numerical 
and/or analytical simulations both of NMR spectra in general and of the NMR emulation 



of quantum computing [|15 



Once proven that the kinematics and the dynamics of the physical systems used in NMR 
spectroscopy are ordinary classical ones, the issue of if and how these systems attain quantum 
computation capabilities is solved. Even though they turn out to behave as classical analogue 
computers, they share with quantum computers, as shown in the following, what up to now 
was considered their unique feature of the exponential growth of state space dimensions with 



the physical dimensions |T^. This answers the question about what makes NMR more 
powerful than conventional classical computing if it does not involve entanglement. 

To be specific the Hamiltonian operator of N identical molecules, each of them containing 
just 2 magnetic nuclei (to keep algebra to the minimum), for times shorter than decoherence 
times, is given by the sum of the Zeeman and the J-coupling terms: 



N 2 N N 

= ^ = - ^ 7p ^ 5 ■ - J ^ ■ (1) 

a=l p=l a=l a=l 

where 5*"^ denotes the spin relating to the pth nucleus in the ath molecule and i? is a 
uniform external generally time dependent magnetic field, by which the Hamiltonian is 
invariant under molecule permutations. Since no external action on the system is going to 
act on different molecules in different ways, one can naturally consider as physical algebra 
the algebra of permutation invariant operators. While the Zeeman term can be written in 
terms of the macroscopic spin vectors = J2a=i "S^"^) this is not so for the J-coupling term, 

by which the tensor operators 

N 

Cjk = E • sf (2) 

a=l 

are to be treated as independent ones, while their commutation relations both between 
themselves and with the macroscopic spin vectors may be computed in terms of their explicit 
expression (Q). To avoid unnecessary complications in illustrating the main ideas, it is 
appropriate to limit consideration to spin 1/2 nuclei, in which case the commutators, after 
some algebra, for ^ = 1 are given by 

= iejkiCii, (3) 

where both commutation relations for spin in general and anticommutation relations for the 
spin 1/2 Clifford algebra 

= iSa/sSpgEiklSf, [S'i^^,S'^^]_^^ = -Sik (4) 

have been used, while of course the macroscopic spin vectors fulfill the same commutation 
relations as the microscopic ones. Here and henceforth summation is implied on repeated 
tensor indices. 

By using these equations and the expression for the Hamiltonian operator as an element 
of the permutation invariant algebra 
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Qap q/Sq 



H = -jZlpB ■ - JCu, (5) 
p=i 



(6) 



one obtains the evolution equations in the physical algebra: 

dCji/dt = i [H, Cji] = —'^2^ikiBkCji — "yiejkiBkCn — Je^i {Sf — Sf) 
dSf/dt = I [H, Sf] = —fpe,,kB,Sj: + i-iyje,,kCk„ p=l,2. 

These equations, as announced, are linear, by which they hold for expectation values 
too. To be specific they correspond to Hamiltonian classical equations with a classical 
Hamiltonian having the same expression as the quantum one in eq. and fundamental 
Poisson brackets obtained by the fundamental commutators in eq. (0) according to the 
reverse Dirac prescription A,B ih{A,B}. Of course the same equations hold for the 
vector and the tensor operators of a s.m., which in turn are algebraically equivalent to the 
Heisenberg equations for a s.m., this showing why the bulk in consideration can emulate a 
s.m. dynamics, and vice versa why the s.m. model works in NMR spectroscopy. 

It should be stressed that this result is independent from the simplifying assumptions 
concerning the single spin value and the number of magnetic nuclei per molecule, even 
though the corresponding equations are in general more complex. Consider for simplicity the 
generalization relevant to most NMR computing experiments, removing only the restriction 
on the number of magnetic nuclei. In such a case the commutators between the tensor 
operators CJ^ = Y.a=i S°'^S'j^'^ give rise to terms linear in third rank tensor operators CJ^[ = 
X]a=i whose evolution equations have to be added to eq.s @. In order to do 

that commutators between third an second rank tensor operators have to be evaluated and 
so on. The end result amounts to consider all tensor operators 

N 

Cf = 5'f , CJ^, C^jli , Cj^jj'i^j^ = X! •••'^in"' ^ — ^I'^i 3j PiQi"^ = 1; 2, ...n, (7) 

a=l 

which generate a (4" — l)-dimensional Lie algebra, since they are linearly independent and 
all their mutual commutators turn out to be linear in terms of them. While the Hamiltonian 
operator retains its structure 
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H = -tl,B.S^-Wj,,C!^, (8) 

p=l p=^g 

the generalizations of eq.s (P), though staying hnear, grow in their number and become 
more complex, except for those giving the time derivatives of macroscopic spin vectors, 
which essentially keep their form. 

As to the initial conditions for the classical evolution equations with reference to NMR 
spectroscopy, they correspond to the expectation values in the thermal state, which can be 



equivalently considered either as a state of the physical algebra ||T^ or of the total algebra 
generated by all microscopic spin operators. In the latter case for the generic tensor operator 
C one gets 

Tr[(E^^,C")exp {-PELiH")] ^ ^ Tr exp (-/jff")] 
Trexp(-/5E^=ii^") k Tr [exp ' 

which obviously coincides with times the expectation value of the corresponding s.m. 
observable with respect to the s.m. thermal state used in the traditional approach. 
In conclusion some remarks are in order. 

The plethora of classical variables at first sight may look puzzling, since most of them are 
not easily accessible. While static measurements only probe, in principle, the three compo- 
nents of magnetization, the very precise dynamic methods of modern NMR spectroscopy 
can give an indirect access, according to the values of the chemical shifts and the J-couplings, 
to all or some of them, which then play the role of more or less hidden variables. Paren- 
thetically, every possible symmetry of the s.m. Hamiltonian, due to the specific values of 
chemical shifts and J-couplings adds together with the permutational symmetry in reducing 
the physical algebra and then the number of classical variables. On the other hand one can 
solve the evolution equations for all variables but for macroscopic spin vectors in terms of 
these ones, whose equations then become nonlocal in time and completely describe the sys- 
tem. Of course the seemingly infinite memory time of these hysteretic equations disappears 
if one takes account of the dissipative terms to be added to describe the long time behavior 

El. 



One of the sources of what I call the mythology of NMR spectroscopy is its explicit 
demonstration of spinor character for a spin- 1/2 nucleus |]12|. While this is considered 
usually as a proof of some sort of quantumness implied in NMR measurements [|1|, the 
present approach shows unambiguously that quantum mechanics is only involved, as for 
most computations in condensed matter physics, in deriving equations connecting classical 
observables. If one considers for instance spin-1 magnetic nuclei, then the operator algebra 
of the single-spin state space is not linearly generated by the spin operators, which leads 
to the appearance of more general tensors involving eight linearly independent operators 
generating (with the identity) the single spin associative algebra Then the value of the 
spin is explicitly reflected in the number of independent classical variables and in the form 
of classical equations. On the other hand this has to be so if the equivalence between the 
differential equation for the s.m. density matrix and the classical bulk evolution equations 
holds true, by which for instance for spin 1/2 just 4" — 1, namely the number of independent 
real elements in the density matrix, classical equations are needed, while, if for instance only 
vector nuclei are present, this number becomes 9" — 1. 

While what has been shown proves the substantial equivalence between quantum com- 
putation and classical bulk computation in systems like the NMR ones, which exhibit an 
exponential growth of the number of classical macroscopic degrees of freedom with the 
number of quantum microscopic ones, this should not raise unjustified optimism about the 
perspectives of this kind of classical non conventional computation. The author has the 
marked feeling that the decrease of the signal-noise ratio as the classical macroscopic de- 
grees of freedom grow in number will prove to be a general obstruction on the way of using 
it as an alternative to quantum computation. Within the present approach the analysis of 
this issue, namely of the signal strength, amounts to the evaluation of the aforementioned 
dissipative terms p5| . 
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